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ABSTRACT 

This paper explores the effect of the LMC on the mass estimates obtained from the timing 
argument. We show that accounting for the presence of the LMC systematically lowers the 
Local Group mass (Mlg) derived from the relative motion of the Milky Way-Andromeda 
pair. Motivated by this result we apply a Bayesian technique devised by Penarrubia et al. 
(2014) to simultaneously fit (i) distances and velocities of galaxies within 3 Mpc and (ii) 
the relative motion between the Milky Way and Andromeda derived from HST observa¬ 
tions, with the LMC mass (Mlmc) as a free parameter. Our analysis returns a Local Group 
mass Mlg = 2.64;|^'3| x IO^^Mq at a 68% confidence level. The masses of the Milky Way, 
TfMw = 1-04;!;q23 X IO'^Mq, and Andromeda, Mm 3 i = x lO'^M©, are consistent with 

previous estimates that neglect the impact of the LMC on the observed Hubble flow. We find 
a (total) LMC mass Mlmc = 0.25^® Qg x IO'^Mq, which is indicative of an extended dark mat¬ 
ter halo and supports the scenario where this galaxy is just past its first pericentric approach. 
Consequently, these results suggest that the LMC may induce significant perturbations on the 
Galactic potential. 

Key words: Galaxy: kinematics and dynamics; galaxies: evolution. 


1 INTRODUCTION 

Despite being the two most luminous satellites of the Milky Way 
(MW hereinafter) the total masses of the Magellanic Clouds re¬ 
main largely unconstrained. One of the key uncertainties is their 
accretion time. Given that satellite galaxies orbiting around a mas¬ 
sive host lose a large fraction of their dark matter halo enve¬ 
lope to tides (e.g. Penarrubia et al. 2008, 2009), an early accre¬ 
tion scenario would truncate the halo profile and limit the to¬ 
tal mass to the mass enclosed within their luminous radius, i.e. 
M ~ M(< /-max) = rmaxvLx/G, where v^ax is the peak velocity of 
the circular velocity curve and rmax the outer-most radius with 
known members. With a circular velocity curve that peaks at Vmax — 
90kms“* and remains flat out to 8.7kpc (van der Marel & Kalli- 
vayalil 2014) the Large Magellanic Cloud (LMC hereinafter) has a 
(minimum) bound mass Mlmc(< 8.7kpc) = 1.7 x lO’^M©. For the 
Small Magellanic Cloud (SMC hereinafter) Vmax — 60kms“' and 
rmax — 3.5kpc (Stanimirovic et al. 2004), which yields Msmc(< 
3.5 kpc) = 2.9 x 10*M©. Taking these estimates at face value and 
adopting the MW mass derived from the dynamics of galaxies 
in and around the Local Group (Mmw = 0.8 x lO'^M©; Diaz et 
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al. 2014; Penarrubia et al. 2014; see also §3) yields a mass ra¬ 
tio (Mlmc-I-Msmc)/Mg — 0.02, which suggests that the combined 
masses of the Magellanic Clouds do not contribute significantly to 
the mass budget of our Galaxy, Mg « Mmw -I-Mlmc -hMsmc- 

This result changes dramatically if the Magallanic Clouds are 
currently past their first pericentric approach to the MW, as sug¬ 
gested by current HST data (Besla et al. 2007; Kallivayalil et al. 
2013). In the first-passage scenario the dark matter haloes of LMC 
and SMC may extend well beyond their luminous radius, such 
that M 3> M(< rmax). How far beyond is currently an open issue 
given that the outer density profile and extent of dark matter haloes 
remain poorly constrained both theoretically and observationally. 
Lacking additional data one may invoke abundance-matching ar¬ 
guments, which suggest that galaxies with the luminosity of LMC 
(M* ~ 2.9 x IO’M©; van der Marel et al. 2002) reside in haloes 
with a virial mass ~ 2 x 1O’*M0 (Mosteret al. 2013). Interestingly, 
A-body simulations of Besla et al. (2012) illustrate the SMC needs 
to complete 2-3 orbits about the LMC within the past ~ 3Gyr in 
order to produce the main components of the Magellanic System 
which, given the high relative speed of the SMC with respect to the 
LMC (~ 130kms“'; Kallivayalil 2013), requires Mlmc > 10'*M©. 

Putting these estimates together and taking into account that 
the present-day Galactocentric distance of the LMC is ~ 50 kpc 
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(McConnachie 2012) yields a striking scenario where the LMC 
mass may be comparable to the mass of the MW at the LMC’s loca¬ 
tion, Mmw(< 50kpc) « 3 X 1 O**M 0 (Xue et al. 2008), and a virial 
radius riMc.vir ~ 1 lOkpc that is approximately twice the distance to 
the MW centre. Such a massive LMC would strongly perturb the 
Galactic potential and affect the constraints on the LG mass derived 
from the timing argument (Gomez et al. 2015). 

This Letter aims to constrain the total mass of the LMC by 
modelling the measured distances and velocities of galaxies in the 
Local Volume (hereinafter LV) as well as the relative motion be¬ 
tween our Galaxy and M31 derived from HST observations (§2). 
Our analysis follows up on the work of Penarrubia et al. (2014), 
who use the equations of motion that govern the dynamics of 
nearby galaxies in an expanding ACDM Universe to model how 
the Local Group (hereinafter LG) perturbs the local Hubble flow 
(§3). Model parameters are fitted using a Bayesian methodology 
outlined in §4. Joint constraints on the parameters of interest are 
given in §5 and the results discussed in § 6 . 


2 THE LOCAL GROUP REST FRAME 

We choose a coordinate system whose origin lies at the LG 
barycentre. Following Diaz et al.(2014, hereinafter D14) and Penar¬ 
rubia et al. (2014, hereinafter P14) the barycentre is found by 
demanding that the total momentum of our Galaxy (G) and An¬ 
dromeda (M31) balances to zero. For simplicity we adopt the same 
vector notation as in D14, where Va-s-b denotes the velocity vector 
of A measured from the reference frame S, such that va->b = -\b^a 
and va^b = va^c +vc^b- 

Below we study the impact of the LMC on the (inferred) mo¬ 
tion of the sun with respect to the LG barycentre. As a first step, let 
us calculate the total momentum of the LG as 

0 = MGVG-»LG+AfM3lVM31-i-LG (1) 

= Mlmc vlmc-»lg + Mmw VmW-»LG + Mm3 1 Vm3 1 -»LG • 

For better comparison with previous work we introduce the param¬ 
eters /„ = Mg/Mm 3 i and fc = Mlmc/Mmw, where Mq ~ Mlmc + 
Mmw- Using Equation (1) it is relatively straightforward to show 
that the velocity vector of the sun with respect to the LG barycentre 
(the so-called solar apex) can be written as 

VM31-»MW , fmfc 

V0-»LG =-r—T-+V0-»MW- ,, r w. , ^ ' VlMC^MW, (2) 

where V0_>mw = (Uq,Vo + Vq,Wq) denotes the velocity vector 
of the sun with respect to the MW centre. We adopt a circu¬ 
lar velocity of the MW at the solar radius Vo = 239 ±5 km s“' 
(McMillan 2011), a Local Standard of Rest {Uq,Vq,Wq) = 
(ILL 12.2,7.2)kms“* (Schbnrich 2012) and a left-hand coordi¬ 
nate system r 0 _>MW = (-8,0,0)kpc. The velocity vectors of M31 
and the LMC with respect to the MW centre are vm 31 ->mw = 
VM31-S-0 + V 0 ^MW and vlmc^mw = vlmc ->0 + V 0 ^mw, respec¬ 
tively. Heliocentric velocity vectors are derived from line-of-sight 
velocities and proper motions. For Andromeda (vios,Vw,Vw) = 
(-300±4,-119.2±36,-6L9±34)kms“' (Sohnetal. 2012; D14), 
LMC (262.2± 3.4,-465 ±5,56± ll)kms^‘ and SMC (145.6± 
0.6,-236± 19,-341 ± 19)kms^' (McConnachie 2012; Kallivay- 
alil et al. 2013). The heliocentric distances of M31, the LMC 
and SMC are rM 3 i^O = 783 ± 25kpc, rLMC ->0 = 51 ± 2kpc and 
rsMC^O = 64±4kpc, respectively (McConnachie 2012). 

Unfortunately, we do not have access to the transverse motion 
of galaxies in the LV owing to their large heliocentric distances. 


Following Karachentsev & Makarov (1996) one can derive a ‘cor¬ 
rected’ radial velocity in the LG rest frame by adding the projection 
of the solar apex onto the heliocentric radial velocity 

V = +V0-»LG • fs-!.0; (3) 

where denotes heliocentric line-of-sight velocities, and 

rj _>0 = rj_> 0 /|r^_> 0 | is a unit vector. At large distances Equa¬ 
tion (3) approximately corresponds to the radial component of the 
velocity vector, i.e. v « Vj_j.lg • f j-»lg. 

Comparison between Equations (2) and (3) shows that the 
presence of the LMC changes the barycentric velocity of nearby 
galaxies by Av^ = -/m/c(l +/m)“‘(l +/c)“'(vlmc^mw • fg^0). This 
quantity can be positive or negative depending on the location 
of the galaxy with respect to the axis that joins the LMC with 
the MW centre, hence a wrong choice for fc will increase the 
scatter in the distance-velocity relation of galaxies randomly dis¬ 
tributed on the sky. Even if the LMC contributes to a small frac¬ 
tion of the Galaxy mass {fc <C 1, /m ~ 1), its large Galactocentric 
velocity, vlmc->mw — 327kms“*, leads to a non-negligible range 
Avc ~ ±/cl60kms“*. Accordingly, the maximum-likelihood value 
of fc is found in §4 by minimizing the scatter in the observed Hub¬ 
ble flow (see Appendix B of P14 for illustration). 

It is straightforward to show that the mass of the SMC 
can be safely neglected in our analysis. Application of the 
Tully-Eisher (1977) relation to the peak velocities of the SMC 
and LMC measured by Stanimirovic et al. (2004) and van der 
Marel & Kallivayalil (2014), respectively, yields Msmc/Mlmc ~ 
(60kms“*/90kms“*)‘' ~ 0.2. Using their measured separation, 
|rsMC-s.LMc| « 24kpc, and relative velocity, |vsmc^lmc| « 
120 kms^*, to balance the clouds momentum leads to a barycen¬ 
tre displaced by « 4kpc and « 20 km s“' from the LMC, which can 
be neglected when compared with its current Galactocentric dis¬ 
tance, rLMC->MW — 50kpc, and velocity, vlmc->mw — 327kms“*. 
Using similar arguments van der Marel et al. (2012b) conclude that 
neither M33 nor M32 change the barycentre of M31 appreciably. 


3 TIMING ARGUMENT & THE LOCAL HUBBLE FLOW 

In an expanding Universe the equations of motion that govern the 
dynamics of galaxies in the outskirts of the LG can be written as 
(Lynden-Bell 1981, Sandage 1986, Partridge et al. 2013, P14) 

r = -^r-H//o^^Ar, (4) 

where r = |rj_>LG| is the distance of a galaxy to the LG barycen¬ 
tre, M = Mg+Mm 31 is the LG mass, Ho = 67 ± 3kms“*Mpc“' is 
the Hubble constant (Efstathiou 2013 and references therein) and 
Ga = 0.686 ± 0.020 is the dark energy density (Planck Collabora¬ 
tion 2013), which yields a Universe age to « 13.7Gyr. 

For galaxies moving on radial orbits the solution to Equa¬ 
tion (4) is univocally determined by measuring the barycentric dis¬ 
tance (r) and velocity (v) at f = to- The isochrone v = v{r,M,to) de¬ 
fined by Equation (4) is known as the ‘perturbed Hubble flow’. In 
this work isochrones are calculated under the assumption that (i) 
the LG potential is dominated by a static MW-M31 monopole, (ii) 
nearby galaxies can be treated as massless tracers and (iii) the LV is 
isolated. Tests in P14 show that both the quadrupole and the time- 
dependence of the LG potential have a negligible impact on the 
observed Hubble flow, while in §4 we use statistical methods to ac¬ 
count for tangential motions that may arise from the self-gravity of 
the tracer population and the presence of galaxy associations in the 
vicinity of the LG. 
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Neglecting the energy term in Equation (4) and adopting 
dM/dr = 0 yields a Keplerian solution to Equation (4) 


r = a(l-ecosJ 7 ) (5) 

' = {^ m ) 


where rj is the ‘eccentric anomaly’, a = l}/[GM{l-e^)\ is the 
semi-major axis, e = l + [2El}/{GMf'} is the orbital eccentricity, 
E=l/ 2v^ - GM / r and L = r x v are the specific energy and angular 
momentum, respectively. We consider orbits where rj = 0 al t = 0; 
?7 = TT at tum-around radius, rape = a(l+e); and 77 = Itt at pericentre, 
Tperi = a(l -e). The radial and tangential velocity components are 


/ fl \esin77 
\GM ) 1 —ecosTj 

/ a \ Vl - A 
\GM ) 1 —ecosTj 


( 6 ) 


Equations (5) and ( 6 ) can be used to measure the combined 
mass of our Galaxy and Andromeda, M = Mmh +Ma, the so- 
called ‘timing argument’ (e.g. van der Marel & Guhathakurta 
2008 and references therein). Position and velocity vectors now 
describe the relative motion between M31 and our Galaxy, i.e. 
r= |rM31->G|, Vrad = Vm 31 -»-G ' fM31-!.G and Vtan = |vM 31 ^G X fM31->G|, 

where fM 3 i^G is a unit vector pointing toward M31 from the 
Galaxy barycentre. Eollowing similar steps as in §2 we find 


VM 31 -»G = VM 31 ^MW- 


fc 

i+fc 


vlmc-»mw, 


(7) 


and similarly for rM3i->G- 

Eig. 1 shows that the LG mass (upper panel) and the orbital 
parameters of the Galaxy-M31 pair (bottom panel) depend on the 
unknown value of fc = Muac/Mm'n- Models that neglect the LMC 
mass, fc « 0, recover the main results found by van der Marel & 
Guhathakurta (2012a), that is a LG mass, M ~ 4.5 x lO'^M©, and 
an orbital pericentre rpai < 20kpc, which imply that our Galaxy 
and M31 will suffer a head-on encounter ~ 4Gyr from now (van 
der Marel et al. 2012b). As discussed in P14, the mass estimate 
increases up to M ~ 5.2 x IO'^Mq if one adopts Vo = 220km s“', in 
agreement with Li & White (2008). 

These estimates change significantly if the mass of the LMC 
exceeds the mass ratio derived from its internal kinematics. In par¬ 
ticular, for fc > 0.2 the LG mass drops down to M ~ 2.8 x lO'^M©, 
in better agreement with the values independently derived by D14, 
M = (2.5 ± 0.4) X IO'^Mq, and P14, M = (2.3 ± 0.7) x lO'^M©, 
while the pericentric distance increases up to rped ~ 100 kpc, indi¬ 
cating a fly-by encounter rather than a head-on collision. Excluding 
models where the apocentric distance (i.e the tum-around radius) is 
smaller than the current separation between our Galaxy and M3L 
i.e. Tapo > 800kpc, imposes an upper bound fc < 0.5. It is thus clear 
that a robust measurement of the LG mass using the timing argu¬ 
ment requires a tight constraint on the (total) LMC mass. Below we 
attempt to measure fc by modelling the dynamics of LV galaxies. 


4 BAYESIAN INFERENCE 

The main goal of this paper is to constrain the individual masses 
of the LMC, the MW and M31 from the measured locations and 
velocities of nearby galaxies taking into account statistical uncer¬ 
tainties in the circular velocity of the sun around the MW centre, 
Vq, as well as in the cosmological parameters Ho and Da- To this 



0 —'—'—'—'—^^^^^—'—'—'—^—'—'—'—'— 

0 0.1 0.3 0.3 0.4 


L ^LMc/^MW 


Figure 1. LG mass (upper panel), and peri- and apo-centre distances be¬ 
tween our Galaxy and M31 (lower panel) derived from the timing argument. 
Equations (5) and (6), as a function of the mass ratio fc = IUlmc/iUmw 
for two values of the MW’s circular velocity at the solar radius (Vq). Note 
that the LG mass M = Mm 3 i +Mq and the impact parameter of the future 
Galaxy-M31 encounter strongly depend on the value of f. 


end we use the following likelihood function (see P14 for details) 

(8) 


TT ^ 

iVi-Vk,,?' 


2(j2 _ 


where S = iM,fm,fc,Vo,h,Q.A,(Jm) is a vector that encompasses 
the model parameters; D and Vh are heliocentric distances and ve¬ 
locities, respectively, and (l,b) Galactocentric coordinates. For a 
given set of 5 we use Equation (4) to find the LG-centric velocity 
V = v(r,M,to), which is then corrected by the apex motion (2) as 
Vi = v(r,,M,ro)-V 0 _>LG • fi -»0 and inserted into the likelihood ( 8 ). 

Equation (4) is solved upon the assumption that galaxies in 
the LV move on radial orbits. Peculiar motions about the Hubble 
flow are accounted for by the hyperparameter Om, which we in¬ 
troduce in the two-dimensional variance of the i-th measurement 
as of = €y i + ej) j(dV/dD)fD.j. i,.^ + a^,. Here ey and cd denote, re¬ 
spectively, velocity and distance errors (see Ma et al. 2013), while 
dV/dD corresponds to the velocity gradient at a given location 
(D,l,b) returned by our model. Tests carried by P14 show that 
marginalizing over a,,, yields unbiased joint bounds on the param¬ 
eters of interest insofar as the tangential velocities that give rise to 
the peculiar motions are randomly oriented on the sky. 

We use a Bayesian inference code, MultiNest (Feroz & 
Hobson 2008), to produce posterior samplings for our parameters 
and measure the error of the evidence. Flat priors are adopted for 
the parameters of interest (M, /„,, fc, and g,„) and Gaussian priors 
for parameters taken from the literature {Vo, Ho and Ha). 

The data consist of a catalogue Ng = 35 systems compiled by 
P14 within barycentric distances 0.78 < r/Mpc < 3.0 and known 
line-of-sight velocities. Most galaxies at r < 0.78Mpc are satellites 
of either our Galaxy or M31 and move on orbits that are poorly 
described by Equation (4). Beyond 3 Mpc three prominent asso¬ 
ciations introduce significant perturbations in the Hubble flow: IC 

















4 Jorge Penarrubia et al. 



Figure 2. Sampling of posterior distributions of the LG mass, M, the mass ratios between the Galaxy and M31, /,„ = Mg/A/m31 and the LMC and the MW, 
fc = Mlmc/^mw- These distributions marginalize over Hq, Ha, Vq and a,„ (see text). Left panels use the likelihood (8) to fit the dynamics of galaxies with 
0.78 < r/Mpc < 3.0. Right panels apply a weighted joint likelihood (9) to fit the local Hubble flow and the motion of the Galaxy-M31 pair simultaneously. 


342/Maffei-I, M81 and Centaurus A/M83. To minimize such per¬ 
turbations P14 only include galaxies whose distance to any of those 
three major associations is larger than 1 Mpc. 

In this work we also construct models that fit simultaneously 
the dynamics of galaxies in the LV and the relative motion of our 
Galaxy and M31. Following Lahav et al. (2000) we consider a joint 
likelihood weighted by the number of objects in each sample, (Ng + 
1) In £ = Aj In £lv + In £tim, where 


In Vlos; Viv, Vw 



ln[ 27 ra,“] 

2 




andx, = {£), vios, vw, vw} denote heliocentric measurements of M31 
with an error cr,. Xi(S) correspond to the model prediction from 
Equations (5) and (6) expressed in a heliocentric frame. 


5 RESULTS 

Fig. 2 displays posterior distributions for each model parameter as 
returned by MultiNest. Note first that the dynamics of galax¬ 
ies in the LV (left panels) do not provide enough information 
to constrain the total mass of the LMC. Indeed, Fig. 2 shows a 
strong degeneracy between the parameters fc = Muac /Mmw and 
f„ = (Mlmc+Mmw)/4Lm3i. Adopting a LMC mass of fc ~ 0.01 
as in previous works yields models where Andromeda is approxi¬ 
mately twice as massive as our Galaxy (/„ — 0.5), thus recovering 
the main result of D14 and P14. Relaxing this assumption, how¬ 
ever, weakly favours solutions where LMC is a factor ~ 10 more 
massive than suggested by its internal kinematics, i.e. fc ~ 0.2, al¬ 
though models where fc = 0.01 cannot be ruled out with high sig¬ 
nificance. Our analysis also returns a constraint on the LG mass, 
M~2 X IO'^Mq, which appears to be insensitive to the mass ra¬ 
tios fc and fm- Fig. 1 shows that this value is a factor ~ 2.2 lower 
than expected from the timing argument with no Magellanic Clouds 


(~ 4.5 X IO'^Mq), and suggests that a massive LMC may help to 
ease this discrepancy. 

Indeed, fitting the dynamics of nearby galaxies and the relative 
motion between our Galaxy and M31 simultaneously (right pan¬ 
els) provides strong support for a massive LMC. We find Mlmc = 
0.25!o:o8(-o:1?) X IO'^Mq at 68% (95%) confidence, which rules out 
models with Mlmc < 0.1 x IO’^Mq at a 95% confidence level. Such 
a large mass is indicative of an extended dark matter halo surround¬ 
ing LMC and supports the scenario where this galaxy is currently 
undergoing its first pericentric passage around the MW. 

Fig. 3 plots the posterior distributions for the individual 
masses of the LMC, Andromeda and the MW. Comparison be¬ 
tween the upper and lower panels shows that the timing con¬ 
straint shifts the mass of the three galaxies to slighlty higher val¬ 
ues, Mmw = l-04!°;23j!o;445 x 10‘^M©, and Andromeda, Mm31 = 
1.33^0 X lO'^M©, while the LG mass increases to M = 

2.64))0 35(^0 03) X lO’^M©. Interestingly, in all our models M31 ap¬ 
pears more massive than the MW, Mmw/Mm 3 i ~ 0.75, although the 
freedom in the value of fc adds uncertainty to this result. 


6 SUMMARY & DISCUSSION 

HST observations indicate that the LMC is slingshotting from the 
MW at a staggering -H327kms“* (e.g. Kallivayalil et al. 2013). 
This implies that even a relatively light LMC may have a mea¬ 
surable impact on the Galaxy barycentre and the masses inferred 
from the timing argument. Using the timing equations we simul¬ 
taneously measure the (total) masses of the LMC, the MW and 
M31 by fitting the locations and radial velocities of a catalogue 
of 35 nearby galaxies (< 3 Mpc) compiled by P14. An upper limit 
Mlmc < 0.31 x lO'^M© is found with 95% confidence. The LMC 
mass is further constrained from below by fitting the dynamics 
of nearby galaxies and the relative motion between the MW and 
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Andromeda, which yields Mlmc = x IO'^Mq at a 

68% (95%) confidence levels. Our method also puts hounds on 
the masses of the MW, Mmw = 1-04^0 23(^0^) x IO'^Mq, and An¬ 
dromeda, Mm3i = 1 X lO'^M©, which are consistent at 
one-sigma level with previous measurements that neglect the im¬ 
pact of the LMC on the observed Huhble flow (D14, P14). Our 
results indicate that LMC represents ~ 1 /5th of the Galaxy mass, 
Mlmc/Mg = (0.19 ±0.05) x IO'^Mq. That is ~ 10 times higher 
than the mass enclosed within its nominal luminous radius, which 
suggests the existence of an extended dark matter halo surrounding 
this galaxy and supports the scenario where the LMC is currently 
undergoing its first pericentric interaction with the MW. 

In addition to the tests of P14 we have inspected two of our 
main assumptions: (i) galaxies can be modelled as point masses 
and (ii) Large Scale Structures (LSS) do not perturb the local Hub¬ 
ble flow (< 3Mpc). Using live A-body models of Gomez et al. (not 
shown here) we find that tidal stripping and dynamical friction only 
affect our timing mass estimates at a < 5% level. This result shows 
little sensitivity to the choice of density profile for the LMC. Also, 
motivated by Libeskind et al. (2015) we have searched for signa¬ 
tures of a velocity shear in our galaxy sample arising from the mo¬ 
tion of the LG through the cosmic web, but found no statistically- 
meaningful deviation between the velocities predicted by the point- 
mass model and the location of galaxies on the sky. The accurate 
description of the local Hubble flow by the timing argument sug¬ 
gests that neglecting the impact of the LSS on the dynamics of 
nearby galaxies does not unduly affect our constraints. 

The results presented here call for theoretical and observa¬ 
tional follow-up work. For example, live A-body experiments are 
needed to understand the perturbations that a massive LMC may 
induce on the Galaxy constituents (disc, bulge and halo). In addi¬ 
tion, extending the work of Li & White (2008) and Gonzalez et al. 
(2014) to the Local Volume will help to further test our method and 
calibrate the dynamical (e.g. virial) masses of the LMC, MW and 
M31 from the measured timing masses. Future kinematic surveys 
of the stellar populations discovered beyond the nominal luminous 
radius of the LMC (e.g. Majewski et al. 2009) will complement the 
mass estimate derived here using cosmological arguments. 
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